Vortex density fluctuations in quantum turbulence 
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Turbulence in the low temperature phase of liquid helium is a complex state in which a viscous 
normal fluid interacts with an inviscid superfluid. In the former vorticity consists of eddies of all 
sizes and strengths; in the latter vorticity is constrained to quantized vortex lines. We compute 
the frequency spectrum of superfluid vortex density fluctuations and obtain the same f~^^^ scaling 
which has been recently observed. We show that the scaling can be interpreted in terms of the 
spectrum of reconnecting material lines. To perform this calculation we have developed a vortex 
tree algorithm which considerably speeds up the evaluation of Biot-Savart integrals. 

PACS numbers: 67.25.dk Vortices and turbulence in superfluid *He 

47.32.C- Vortex dynamics 

47.27.Gs Homogeneous isotropic turbulence 
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Current theoretical and experimental work explores the relation between turbulence in an ordinary (classical) fluid 
and turbulence in the quantum phases of ^He, '^He and atomic Bose-Einstein condensates. Quantum turbulence[l| 
shares important features with classical homogeneous isotropic turbulence the most important is the Kolmogorov 
^-5/3 g]2g]-gy spectrum[3| where k is the wavenumber. 

^He consists of two components: an inviscid superfluid components (associated to the quantum ground state) and 
thermal excitations, which make up a viscous normal fluid component. What makes helium particularly interesting 
is that superfluid vorticity is concentrated in line singularities of flxed circulation k = h/m, where h is Planck's 
constant and m is the mass of one helium atom (in '^He the relevant boson is a Cooper pair). Normal fluid vorticity is 
unconstrained, as in classical flows. Turbulence in *Hc is thus a complex doubly-turbulent regime, in which a viscous 
fluid interacts with discrete inviscid vortex lines. 

The intensity of quantum turbulence is characterized by the vortex line density L (vortex length per unit volume) . 
In a striking experiment, Roche et al.^ measured the fluctuations of L in turbulent ^He. They observed that the 
frequency spectrum scales as f~^^^, where / is the frequency. The same scaling was observed in turbulent ^He[^. The 
rapid decrease of the spectrum is surprising because, if one interprets L as a measure of the rms superfuid vorticity 
{ujs = nL), it seems to contradict the classical scaling of vorticity [7| expected from the Kolmogorov energy spectrum, 
which increases with k. The aim of this letter is to shed light onto this problem. 

Since the typical distance between superfluid vortices I w L~^l'^ is orders of magnitude bigger than the vortex core 
radius ao, we model vortices as space curves s(^,<) where ^ is arclength and t is time. The curves are numerically 
discretised by a large, variable number of points (j = 1, • • • N^. The governing equation is Q 



— = Vs + as' X (v„ - Vs) - a's' x [s' x (v„ - v^)] , (1) 

where s' = dsjd^ is the unit tangent vector at s, a, a' are temperature dependent friction coefHcients @, v„ is the 
normal fluid's velocity, and the velocity Vs which the vortex lines induce on each other at the point s is given by the 
Biot-Savart (BS) law: 



The line integral extends over the entire vortex configuration H. Vortex lines reconnect I lOi when they become 
sufficiently close to each other, provided that the total length (as a proxy for energy) is reduced[ll|. The discretization 
technique is standardjT^: the reconnection technique and the de-singularization of the BS integral are described 
elsewhere [Tsj. 

The difficulty of this vortex filament method is the computational cost of the BS law which scales as N"^ (the 
velocity at one point depends on an integral over all other iV — 1 points) ; this prevents calculations of intense vortex 
tangles (large N) for sufficiently long times to make realistic comparison with experiments. The same difficulty 
arises in astrophysical N-body simulations (the force of gravity on one body depends on the other A'^ — 1 bodies); in 
this context, the problem was solved by the development of tree algorithms whose computational cost scale as 
7Vlog(A^) with small loss of accuracy [isl. [T6j. 

To achieve our aim we have developed the following tree algorithm for vortex dynamics. At each time step, points 
are grouped in a hierarchy of cubes which is arranged in a three-dimensional oct-tree structure. We construct the 
tree top down, first dividing the computational box (root) into eight cubes, and then continuing to divide each cube 
into eight 'children', until either a cube is empty or only contains one point. As we create the tree, we calculate the 
total vorticity contained within each cube and the 'center of vorticity' of the cube from the points that it contains. 
The time required for constructing the tree scales as iVlog(A^), so it feasible to 'redraw' the tree at each time step. 
Fig. mniustrates this procedure in two dimensions. 

To calculate the induced velocity at each point s we must 'walk' the tree, and decide if a cube is sufficiently far. 
This is done using the concept of the opening angle Q [ij] (as corrected by Barnes[l^[l3 to avoid errors if the center 
of vorticity is near the edge of the cube). Let w be the width of the cube, d the distance of the center of vorticity from 
s and (5 the distance from the center of vorticity and the geometrical center of the cube. If(i>w/0 + (5we accept 
the cube, and its contribution is used in computing the velocity via Eq. [51 If not, then we open the cube (assuming 
it contains more that one point) and repeat the test on each of the child cubes that it contains. The tree-walk ends 
when the contributions of all cubes have been evaluated. 

We tested the tree algorithm up to iV = 2000 points (practical limit of the BS law) using different values of 0. We 
verified the N log(iV) scaling for both the construction of the tree and the calculation of the total velocity. We found 
that the relative deviation of the velocity computed via the tree algorithm from the exact BS velocity is at the most 
0.25% \i d — 0.7, which we take as the critical opening angle hereafter. 
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FIG. 1. (Color online) Illustration of the tree construction in two dimensions (quad-tree). The points (red dots) are inclosed 
in the root cell (a), which is divided into four cells of half size (b), until (c,d) there is only one point per cell. 




FIG. 2. Vortex line density L (cm ^) vs time t (s) corresponding to Re„ — 22.7 (dot-dashed line), Re„ — 57.1 (solid line), and 
Ren = 112.9 (dashed line). 



The computational box is a cube of size D = 0.075 cm with periodic boundary conditions. When evaluating the 
BS integral, for each point in the box we consider the other 3'^ — 1 = 26 boxes around it; this periodic wrapping is 
easily obtained using the tree structure.. 

To model the turbulent normal fluid of the experiment^ we use a Kinematic Simulation (KS) [l^, in which the 
normal fluid velocity at position s and time t is prescribed by the following sum of M random Fourier modes: 



M 

i) = ^ (A™ X k„i cos (j>jn + B„i X k™ sin 



(f'ni) , 



(3) 



ni—l 
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FIG. 3. Saturated vortex tangle at t = 2.0 s with N = 55, 359 and L — 91, 733 cm ^ , corresponding to Re„ — 507. 




FIG. 4. (Color online). Normal fluid's (upper two lines) and superfluid's (lower two lines) energy spectrum E{k) vs k (cm ^). 
Grey lines: Re^ = 112.9 (L = 8889 cm'^ ki = 27r/^ = 592 cm"^); black lines: Ren = 49.85 (L = 7058 cm'^ ke = 527 cm'^). 
Dashed lines: k~^^^ (top) and k"^ (bottom) scalings. 
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FIG. 5. (Color online) PDF of superfluid velocity components (cm/s) i}sx (blue circles), Vsy (red triangles) and Vsz (green 
asterisks) sampled over the vortex-points for Re„ = 112.9. The overlapping black dotted, dash dotted and solid lines are 
respectively the Gaussian fits to the same data, gPDF(iisi) = (l/((T\/27r)) exp {—{vsi — fi)'^ /{2a^)), [i = 1, 2, 3) where a and p. 
are the standard deviation and the mean. 
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FIG. 6. (Color online). Power spectral density of fluctuations of L (arbitrary units) vs / (s ^) at t = 10 s corresponding to 
Re„ = 507 as in Fig. [3] The best fit to the data is f^^'"^^. The dashed line shows the f^^^'^ scaling. 




FIG. 7. (Color online). Power spectral density of fluctuations of length of reconnecting material lines (arbitrary units) vs / 
(s^^) sX t — 10s, corresponding to Re„ = 33.28. The best fit to the data is /~^ ''*. The dashed line shows the f^^^^ scaling. 



with = km -s+Wmt, whcrc km and a;„i = y^k^E{km) are wavevcctors and frequencies. Via an appropriate choice of 
and Bm, the energy spectrum of v„ reduces to the Kolmogorov form E{km) oc km^^ for 1 <C A: fc^/, with k = 1 
at the integral scale and fcjvf at the cut-off scale. The effective Reynolds number Re„ = (fc^/fci Y^'^ is defined by the 
condition that the dissipation time equals the eddy turnover time at fc = kM- Like some [l9| previous implementations 
of KS, we have adapted Eq. [3]to periodic boundary conditions. In summary, is solenoidal, time-dependent, and 
satisfies the main properties of homogeneous isotropic turbulence, from the energy spectrum to two-points statistics. 

We use parameters which refer to ^He: circulation n ~ 9.97 x 10^"* cm^/s and vortex core radius oq = 10^^ cm. We 
choose T = 2.164 K (a = 1.21 and a' = —0.3883), larger than in Ref. in order that the back reaction of the vortex 
lines on the normal fluid is negligible and Eq. [3]is justified. Our calculation has thus two independent parameters: T 
and Re„. The initial condition consists of 16 straight vortices at random positions and orientations. 

Figure [5] shows time series of the vortex line density at three different values of Re„. In each case the initial growth 
is followed by saturation to a statistical steady state in which L fluctuates around a mean value. An example of 
a saturated vortex tangle is shown in Fig. |3l By harnessing the power of the tree algorithm wc have performed 
calculations with up to TV = 400, 000 points. 

We construct a 512^ mesh in the xy— plane at the center of the box. At each mesh point we calculate Vg and 
v„, using the tree approximation to the BS integral and Eq. [3] respectively. The corresponding energy spectra for 
two values of Re„ are shown in Fig. ID The fc~^/'^ Kolmogorov spectrum of the normal fluid is clearly visible. The 
superfluid follows the Kolmogorov scaling in the inertial range 1 ^ /c <C fcjvf, in agreement with experiments [3] . In 
the range k > /cm the normal fluid is essentially at rest, and the friction dampens Kelvin waves on quantized vortices 
(preventing a cascade of energy to larger k which would happen if a = a' = 0, as discussed in Ref. [l3|); a k~^ scaling, 
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typical of individual straight vortex lines, is visible in this range. 

Despite the classical nature of the superfluid energy spectrum, the statistics of superfluid velocity components 
display power-law behavior. The probability density functions (normalized histogram, or PDF for short) scale as 
PDF{vs.i) oc {i = 1,2,3) with average exponent b = —3.1, see Fig. [5j This scaling was observed in turbulent 
helium experiments [20j . and was calculated in turbulent atomic condensates (2l]| and helium counterflow j22) : its cause 



is the singular nature of the superfluid vorticity 2lJ . The vortex line velocity ds/dt obeys non-Gaussian scaling too. 
The statistics of velocity components in ordinary turbulence, on the contrary, are Gaussian (23|. 

Finally we compute the frequency spectrum of the fluctuations of the vortex line density in the saturated state. 
Fig. [6] shows that the spectrum scales as f~^^^, for large /, as observed in the experiments [l, Hi- The fact that the 
high frequency regime where this scaling takes place is not exactly the same as in the experiment is less important, 
and arises from our choice of temperature, hence the value of saturated L. We observe the same f~^^^ scaling at all 
values of turbulent intensities Re„. 

What is the reason for this scaling ? Roche et al. 0, [1^ argued that the more randomly oriented vortex lines 
(which particularly contribute to line length and second sound attenuation) have some of the statistical properties of 
passive scalars. To test this idea we perform numerical simulations in which vortex filaments are replaced by passive 
material lines which evolve according to ds/dt ~ v„ (wc do not switch off the rcconnection algorithm, otherwise the 
vortex length would grow indefinitevely) . We find that the length saturates at values larger than the vortex line 
density and that the spectrum of the fluctuations scales again as ,f~^^^, as in Fig. [S] Being L positive definite p6l|. 
there is no conflict with the vorticity spectrum. 

In conclusion, our calculations reproduce the main observed features of quantum turbulence: (i) the classical 
k^^/^ scaling of the energy spectrum observed at large scales by TabelingQ (thought to be associated with large- 
scale, energy-containing polarization of vortex linesQ); (ii) the observation of non classical (non Gaussian) velocity 
statistics [131 (macroscopic manifestation of singular vorticity[2l|); and (iii) the f~^^^ spectrum of the vortex line 
density fluctuations observed at large frequency Q. The natural question is whether they are independent. Our 
results also support Roche's interpretation[23| that vortex density fluctuations arise from random vortex lines which 
behave as reconnecting material lines (while most of the tangle's energy is in the large scale motion). 

The vortex tree algorithm could be further speeded up by parallelizationfi^. Its power should allow us to tackle 
problems which require large A^, for example the detection of anomalous scaling |28{ and, in the T = limit, the 
transition from the Kolmogorov energy cascade at small k to the Kelvin waves cascade at big k [1^ . 
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